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We solve the Klein-Gordon equation for a massive, non-minimally 
coupled scalar field, with a conformal coupling, undergoing cosmo- 
logical evolution from a radiation-dominated phase to a future sud- 
den singularity. We show that, after regularisation, the energy of 
the created particles is zero and the back-reaction from quantum 
effects does not change the evolution of the universe near the future 
singularity and cannot prevent the finite-time sudden singularity. 
PACS number: 98.80-k 



1 Introduction 

The combined data from the anisotropy of the cosmic microwave background 
radiation, supernovae type la, baryonic acoustic oscillations and matter power 
spectrum indicate that the dominant component of the mass-energy content 
of the universe must be gravitationally repulsive and driving the acceleration 
of the universal expansion. In general, such a component can be represented 
by a fluid with an equation of state p — ujp, with to < — i. If a; is constant 
and the spatial sections of the universe is flat, the recent results of 7-years 
WMAP observations indicate that ui = —1.10 ± 0.14[T] at la. This indicates 
that the null energy condition is violated, and the expanding universe may hit a 
universal singularity in its future evolution, after a finite proper time. This type 
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of future singularity has been named " big rip" , when there is a divergence in 
the density of the exotic dark energy fluid, while the scale factor goes to infinity 
[2]. This strange behaviour is connected with the violation of the null energy 
condition. The possibility that the null energy condition (p + p > 0) is violated, 
and cosmological singularities may occur in the future evolution of the universe, 
has been evoked frequently in the literature in the last years. The possibility of 
future singularities was pointed out for the first time in reference [3] and their 
occurrence does not require the violation of the energy conditions. There is 
also a later discussion in ref [7] ■ A recent example of future singularities which 
does not violate the null energy condition is the 'big brake' singularity, which 
emerges from the DBI action [?J [5] . The big brake singularity has the curious 
property that it can be traversed by a pointlike particle [6]. 

A milder type of finite time singularity is the so-called "sudden singularity". 
The sudden future singularity occurs without violation of any energy condition 
(so p + p > and p + 3p > at all times) . This singularity is characterized by 
a finite value for the scale factor, its first time derivative, and for the density, 
while the second derivative of the scale factor and the pressure diverge at finite 
time [3]. They are singularities of the weak sort discussed by Tipler [8] and 
Krolak 0. 

In general, it is believed that the fate of the universe near a singularity 
(past or future) must be affected by quantum effects arising in the extreme 
conditions that exist in its spacetime neighborhood. In the case of the big rip, 
this problem has been treated for example, in references [101 [TTJ [12]. In these 
investigations, it was found that the quantum effects are important. But, the 
conclusions concerning the back-reaction of the quantum effects on the evolution 
of the universe were harder to decide unambiguously. In the case of the sudden 
singularity, quantum effects were studied in references [T31 [H] , and the results 
indicated that quantum effects did not change the evolution near the singularity. 
These results were all obtained for a massless scalar field. 

It is important to extend these results to the case of massive scalar fields. 
Here, we extend earlier studies to the case of a massive scalar field, non- 
minimally coupled to gravity, with a conformal coupling parameter, in a universe 
which possesses a sudden singularity when quantum effects are absent. We will 
show that the singularity is unscathed by the back reaction from the quantum 
field. The complexity of the background evolution of the sudden singularity 
restricts the possibility of a complete analytical solution. Hence, we simplify 
the model by considering two phases: a primordial radiative phase and a sud- 
den singular phase. Under these conditions, we can determine the evolution of 
the quantum fields. In particular, we can calculate the Bogoliubov coefficients 
and determine the energy density of the particles created by quantum effects 
near the singularity. After regularisation, we find a null result implying that 
the sudden singularity persists, unaffected, by the quantum effects due to the 
massive scalar field. 

The paper is organized as follows. In the next section, we determine the 
master equation for the massive, non-minimally coupled, scalar field. In section 
III, we set out the cosmological background containing a sudden singularity, 
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and in section IV, the master equation is solved. The Bogoliubov coefficients 
are determined in section V. In section VI, we regularise the expression for the 
energy density, and show that the final result is zero. In section VII we discuss 
our conclusions. 



2 The master equation 

We shall investigate quantum creation of massive particles near a sudden cosmo- 
logical singularity at finite time. We aim to determine if these quantum effects 
are strong enough to change the evolution of the universe near a classical sudden 
singularity. In order to do this, we must construct a specific model. We will 
concentrate on the evolution of a massive scalar field, <ft, non-minimally coupled 
to gravity, described by the following Lagrangian: 

L = \m 2 <t> 2 + ^R4> 2 - (1) 

The presence of mass takes us one step further than previous work concerning 
particle production near sudden singularity [13] , which considered only massless 
scalar particle production. The non-minimal coupling should strengthen any 
quantum effects near the singularity, and also introduces some technical features 
that facilitate obtaining exact solutions to the problem. 

By variation of the Lagrangian with respect to <fi,~we obtain the field equation 

R 

□0 + to 2 0--0 = O. (2) 

Under variation with respect to the metric, the Lagrangian also gives the momentum- 
energy tensor: 

2 1 1 2 2 4> ( \ 1 2 



where 



is the Einstein tensor. 

If we use equation ([2]) together with the flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) metric for the universe, 

ds 2 =dt 2 -a 2 (t)(dx 2 +dy 2 + dz 2 ), (3) 



where a(t) is the expansion scale factor, we find 



where over-dots denote derivatives with respect to the comoving proper time, 
t. This equation can be rewritten in terms of the conformal time, 77, defined by 
dt = adri, as 

0" + 2 - ( f ) >+( k 2 + ma 2 + —)(/) = 0. 
a \ a J 

The primes denote derivatives with respect to the conformal time. 
Redefining the scalar field as, 




we obtain the following equation for 

X" + (fc 2 +mV)x = 0. (4) 
This is the principal equation we will work with from now on. 



3 Two cosmological eras 

The sudden singularity can be described by the following expression for the 
expansion scale factor (TSJ [16j 



a(t)=(l)V-l) + l-(l-^) B 1 (5) 

where t s is the time where the sudden singularity occurs, and a s is the value of 
the scale factor at this moment. Moreover, < q < 1 and 1 < n < 2 where q 
and n are free constants and no relation is assumed between the pressure p and 
the density p. We have two asymptotic phases: 



• Primordial phase, t 0: 



a ( - 1 (o. - 1), (6) 

/ + \ 9- 2 
1 I 1W * 



i(g-l)fJl) (o.-l). (8) 



• Singular phase, t—}t s 



f(« s -i), (10) 



a -> a s , (9) 
n 

n 1 
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There is a radiation-dominated primordial phase if q = 1/2. On the other 
hand, in the singular phase, the scale factor and its first derivative approach 
constants, and the second derivative, a, diverges as t — > t s , since n < 2. 

In terms of the conformal time, dr/ = a~ 1 dt, we have for the scale factor 
evolution to leading order: 

• Radiation phase: 

a = ciQr). (12) 

• Singular phase: 

a = a s . (13) 

The scale factor and its first derivative must be continuous during the tran- 
sition of one phase to another. If a s is the scale factor value at the moment of 
the transition, and Hq the corresponding Hubble parameter, then the transition 
moment is given by r/t = l/(-Ho a s) and ao = Hq a 2 . 

The isotropic and homogeneous form we have assumed for the cosmological 
evolution of the scale factor, a{t), towards a sudden singularity captures the 
essential features of the general solution near such a singularity. On approach 
to the sudden singularity as t — > t s , the FLRW solution has the linear asymptotic 
form 

a -> a s + q(l - <x s )(l - ^-). 
We can be generalise it to an inhomogeneous metric of the form 

ds 2 = dt 2 - (a a p + Tb aj3 + T n c af 3 + ....) delete' 3 (14) 

on approach to a sudden singularity at r = t — t s = 0, where a Q( g, b a p and c a p 
with a, /3 — 1,2,3 are functions of the space coordinates and n is a constant 
such that 1 < n < 2. In the absence of an equation of state, nine components of 
the symmetric a a p,b a p and c a p tensors are left independent and arbitrary by 
the field equations as r — > and so (fl"4")) is characteristic of part of the general 
solution of the Einstein equations in the vicinity of a sudden singularity [18] . 
A stability analysis of the FLRW solution has also been performed by Barrow 
and Lip [T5]. For other studies of sudden singularities of this type, in general 
relativity and related theories of gravity, see refs [20 1 [2Tj |2"2l 133 1 1^ 1 [23j [2^1 [27 1 

nn usi sni eh isa - 

4 The solutions for the master equations 

Let us return to the master equation First, note that if the mass is zero, 
m = 0, then during both phases the equation becomes, , 

X " + k 2 X = 0. (15) 
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The solution is the same during these two phases, and they can be written in 
the form of plane waves. Since the solution must be continuous, there is no final 
effect in the singular phase, for the trace anomaly or for particle production. 

If the mass is non-zero, the equation during the singular phase takes the 
form, 

X " + (k 2 +m 2 a 2 s ) X = 0. (16) 

The solution is still of plane-wave form, but with a frequency that is affected 
by the presence of the mass term. During the radiative phase, the equation is 

X " + (k 2 +m 2 a 2 r 1 2 ) X = Q. (17) 

which can be re-written as 



where x — ^2maoi], uj = k 2 /2mao, with dots now denoting derivatives with 
respect to x. This is a parabolic cylinder equation, and the solutions can be 
written as, 



X = e 



ci 1F1 r, s, i — +c 2 x 1F1 p, q, i — 



(18) 



where the \F\ are confluent hypergeometric functions (Kummer functions), with 

w 1 1 

r = l 2 + A ' S= 2' (19) 
w 3 3 

p = 1 2 + a ' s= r ( 0) 

where ci,2 are constants. 

The confluent hypergeometric functions may be expressed as a series: 

.fl(r...«)-l + J. + i^^ + .... 

The series represented by the constant c\ is even, while the solution represented 
by the constant c-i is odd. 

Let us re- write the solution as, 

X = ciXi+c 2 X2- (21) 
Considering the series expansion described above we find, 

» - 1 -4-i(n-T)*' + -- < 22 > 
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An important detail is that if we want to take into account the mass in equation 
p7|) we must consider the series at least until order x 5 (x 4 for the even series), 
while if the mass is not taken into account, we can stop at order x 3 (x 2 for the 
even series). The initial conditions are imposed for 77 — > (x — > 0) when the 
mass is negligible. Considering the expansion for the massless case, and going 
back to the notation in terms of the conformal time, we find 

/ k 2 rf\ sj2ma^ ( k 3 T] 3 \ . . 
X = C1X1 + C2X2 ~ ci 1 — +c 2 : kr] — (24) 



~ C1COSK77 + C2 r SinK77. (25) 

K 

Choosing 

Cl = Wk ' C2 = 7^7F (26) 

the solution for 77 — > can be written as, 

X ~ ^Le- 4fc ", (27) 



2k 

corresponding to the vacuum initial state. 

With those choices for the constants ci and C2, we now have the solution at 
any time, also for the massive case: 

X=^=e~ 1 ^ iFi(r,s,i% r ^]+iy/axiF 1 (p,q,i% r ^\ , (28) 



'2k L V 2 / 
where we have restored the mass in the expressions 



5 The Bogoliubov coefficients 

Suppose that we approximate the solution during the radiation era by a mass- 
less field. This is equivalent to assuming that we have the massive solution 
shown above, but in the limit 77 — >• 0. Hence, we have the following solutions 
corresponding to the two phases: 

4>k(r]) = — (primordial phase), (29) 



>2k 

M = eoie^+foae-* 5 " (singular phase), (30) 



where £01,02 are constants and Cj — \Jk 2 + m 2 a^. 

Now we impose the matching conditions for these two solutions at rj = rj c 
by requiring continuity of the function and of its first derivative. We obtain the 
following relations: 

£01 - i-^fl + ^V* -05 " . ( 31 ) 



2 V2k V " 

1 1 / * 

2 V2k V u 



£ 02 - I -±= (l-^) e ^ k+ ^. (32) 
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We notice that when m = 0, £01 = l/v2fc and £02 = and so the solution is 
the same in the two phases and there are no particle production effects. 

On imposing the quantisation, we find the following expressions for the two 
phases: 

<f>k(v) — a H j^ a ^ (primordial phase), (33) 

Ml) = Koie^ + foe-^Ja 

+ (Sie^+^e^ja 1 (singular phase), (34) 

where a and a' are the creation and annihilation operators. These solutions 
(and their derivatives) are continuous at r\ = rj c . They can be re-written as, 



V2k~ ' V2k' 



9k(rf) — ^ a H i^ a ^ (primordial phase), (35) 



<j)k(v) = b + - (singular phase). (36) 

Hence, we have 

b = \/2(i(£oia + £o2a t ). 

In this way we find expressions for the Bogoliubov coefficients that connect the 
quantum modes during the different phases |34j : 



a = V2w£oi , /3 = V2w£o2- 
The normalisation condition, 

aa* - /3(3* = 1, (37) 

is satisfied. When the mass is zero (ie a conformally coupled scalar field), we 
have a = 1 and (3 — 0. The coefficient (3 is associated with the created particles 
[3"4"] . Hence, the number of created particles for each mode k is, 

1/ fc^ 2 



N* = PP = Z[l-~) • ( 38 ) 
while the energy of each mode is 

Pk = kN k . (39) 
An integration over all fc-modes gives, 

fc^ 2 



p = J p k d- i k = Tr J fc 2 w^l--rj dk. (40) 

This expression clearly diverges so it is necessary to regularise it. But, heuristi- 
cally, since it is a polynomial expression, it seems clear that after regularisation 
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we must obtain zero. Hence, the particle production should not contribute to 
the energy-momentum tensor and the sudden singularity is unaffected by these 
quantum effects 

Note that the integral (|4T))) admits an analytical solution: 



Pk d 3 k = TT / fc 2 <i(l - -J dk = vr<j k^/W 
2[k+Vk 2 + '. 



rn 
T 



2 



- 2 

m , 
— In 
4 



(41) 



with m = maQ. There is no infrared divergence, but there is a logarithmic 
divergence when k — > oo (ultraviolet limit). 



6 Regularising the energy 

In order to regularise the expression of the energy, we use the n-wave method 
expounded in the reference [35 . This method is based on the Pauli-Villars 
technique used for quantum field theory in Minkowski space-time. First, let us 
write the energy as, 

/•OO 

p= I p k (k,m)k 2 dk. (42) 
J o 

Let us define, 

pT = —pk(nk,nm), (43) 
n 

where n is a parameter that characterizes the order of the divergence. From 
this expression we construct the quantities, 



n— >oo 



cPp< k n) 
d(n~ 



K = ^i 7 ^-. (44) 



The expression for the regularised energy is given by, 

p r k e9 = p k -E° k -El-lEl (45) 

where E® eliminates the logarithmic divergence, E\ the quadratic divergence, 
and the quartic divergence - all those that are normally present in the 
energy-momentum tensor. This regularisation of the energy corresponds to a 
full renormalisation of the coupling constants, as described in [36} I37j. 
We have, 

Pk = Vk 2 +fh 2 -2k+-^=. (46) 



It follows that 

» _ 



Pi n, =p k - (47) 
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(48) 



As we suspect, the renormaliscd energy is zero. There is no effect, and the 
quantum phenomena associated with the cosmoiogicai dynamics do not change 
the character of the sudden singularity or prevent its occurrence. 

7 Conclusions 

In this work we have investigated the fate of the universe near a future sudden 
singularity due to quantum particle production effects by a massive scalar field 
that is non-minimally coupled to gravity. We have used the method of cal- 
culating Bogoliubov coefficients and the expression obtained for the energy of 
the created particles near the singularity is divergent. This divergence is cured 
using the standard n-wave method and the final result after regularisation is 
exactly zero. Hence, we can conclude that the sudden singularity is robust 
against quantum effects due to the presence of a massive scalar field. We have 
used a simple description for the background cosmoiogicai evolution towards a 
sudden singularity at finite time that shares the same time evolution as part of 
the general solution of the Einstein equations near such a singularity. 

An interesting consequence of this result concerns the trace anomaly. A 
massless scalar field conformally coupled to gravity is a particular case of the 
problem studied here. For this particular case, the result is the same, with no 
quantum effects. Hence, the trace anomaly is absent near the singularity, as 
had been speculated in reference [13] . It must also be stressed, concerning this 
massless limit, that the same regularisation method can be used since there is 
no infrared divergence even when the mass is zero. 

Others fields should be considered to test the robustness of a classical sudden 
singularity with respect to quantum effects. However, since the scalar field and 
its first derivative are both constant at the sudden singularity, the Klein-Gordon 
equation implies that in general there will only be a change in the frequency 
of the quantum modes as the sudden singularity is approached. This frequency 
change alone is unable to create influential quantum effects. Hence, we can 
surmise that perhaps the result found here and in references [131 114] can be 
generalised further, and they can be applied to other types of fields. 

One such a generalisation arises if we consider a non conformal coupling, 
with a Klein- Gordon equation given by 



For the radiative phase, the scenario is the same as that studied above, since 
a" = 0. For the sudden singularity phase, the situation is more involved since 



Ucj) + m 2 <f) - £R(j) = 0. 



(49) 



For a flat FLRW metric, this equation reduces to 




(50) 
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there is a singularity in a". However, a transformation of the type <j) — a~ 6 ^x 
(which reduces to our previous transformation for £ = 6 - the conformal cou- 
pling) can eliminate this singularity, leading to a regular equation with a damped 
(anti-damped) harmonic oscillator equation for £ > | (£ < |). This new term 
seems to be harmless since the dissipation (anti-dissipation) effect last only for 
a finite time until the singularity is reached. Hence, we expect that our previous 
results hold even in this more general case. 

Another way to consider the problem of quantum avoidance of future singu- 
larities is to solve to Wheeler-de Witt equation in universes containing future 
singularities. For big rip and big brake singularities, this question has been 
analysed in references [551132], with some indications that, at least for the big 
brake case, the singularity can be avoided due to quantum effects. For the big 
rip case the situation is less clear, see also references [40l [41]. It is a natural 
step to perform such an analysis for the sudden singularity case, an issue we 
hope to consider in the future. 

Acknowledgements: We thank CNPq (Brasil) for partial financial support. 
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